Discovery of Novel Anti-Resistance AR Antagonists Guided by Funnel Metadynamics Simulation.

Androgen receptor (AR) antagonists are widely used for the treatment of prostate cancer (PCa), but their therapeutic efficacy is usually compromised by the rapid emergence of drug resistance. However, the lack of the detailed interaction between AR and its antagonists poses a major obstacle to the design of novel AR antagonists. Here, funnel metadynamics is employed to elucidate the inherent regulation mechanisms of three AR antagonists (hydroxyflutamide, enzalutamide, and darolutamide) on AR. For the first time it is observed that the binding of antagonists significantly disturbed the C-terminus of AR helix-11, thereby disrupting the specific internal hydrophobic contacts of AR-LBD and correspondingly the communication between AR ligand binding pocket (AR-LBP), activation function 2 (AF2), and binding function 3 (BF3). The subsequent bioassays verified the necessity of the hydrophobic contacts for AR function. Furthermore, it is found that darolutamide, a newly approved AR antagonist capable of fighting almost all reported drug resistant AR mutants, can induce antagonistic binding structure. Subsequently, docking-based virtual screening toward the dominant binding conformation of AR for darolutamide is conducted, and three novel AR antagonists with favorable binding affinity and strong capability to combat drug resistance are identified by in vitro bioassays. This work provides a novel rational strategy for the development of anti-resistant AR antagonists.

On the Dynamic Mechanism of Long-Flexible Fatty Acid Binding to Fatty Acid Binding Protein: Resolving the Long-Standing Debate.

Fatty acids (FAs) are one of the essential energy sources for physiological processes, and they play a vital role in regulating immune and inflammatory responses, promoting cell differentiation and apoptosis, and inhibiting tumor growth. These functions are carried out by FA binding proteins (FABPs) that recognize and transport FAs. Although the crystal structure of the FA-FABPs complex has long been characterized, the mechanism behind FA binding and dissociation from FABP remains unclear. This study employed conventional MD simulations and enhanced sampling technologies to investigate the atomic-scale complexes of heart fatty acid binding proteins and stearic acid (SA). The results revealed two primary pathways for the binding or dissociation of the flexible long-chain ligand, with the orientation of the SA carboxyl head during dissociation determining the chosen path. Conformational changes in the portal region of FABP during the ligand binding/unbinding were found to be trivial, and the overturn of the ″cap″ or the unfolding of the α2 helix was not required. This study resolves the long-standing debate on the binding mechanism of SA with the long-flexible tail to FABP, which significantly improves the understanding of the transport mechanism of FABPs and the development of related therapeutic agents.

Identification of Monobenzone as a Novel Potential Anti-Acute Myeloid Leukaemia Agent That Inhibits RNR and Suppresses Tumour Growth in Mouse Xenograft Model.

Acute myeloid leukaemia (AML) is one of the most common types of haematopoietic malignancy. Ribonucleotide reductase (RNR) is a key enzyme required for DNA synthesis and cell proliferation, and its small subunit RRM2 plays a key role for the enzymatic activity. We predicted monobenzone (MB) as a potential RRM2 target compound based on the crystal structure of RRM2. In vitro, MB inhibited recombinant RNR activity (IC50 = 0.25 µM). Microscale thermophoresis indicated that MB inhibited RNR activity by binding to RRM2. MB inhibited cell proliferation (MTT IC50 = 6-18 µM) and caused dose-dependent DNA synthesis inhibition, cell cycle arrest, and apoptosis in AML cells. The cell cycle arrest was reversed by the addition of deoxyribonucleoside triphosphates precursors, suggesting that RNR was the intracellular target of the compound. Moreover, MB overcame drug resistance to the common AML drugs cytarabine and doxorubicin, and treatment with the combination of MB and the Bcl-2 inhibitor ABT-737 exerted a synergistic inhibitory effect. Finally, the nude mice xenografts study indicated that MB administration produced a significant inhibitory effect on AML growth with relatively weak toxicity. Thus, we propose that MB has the potential as a novel anti-AML therapeutic agent in the future.

A Novel In Situ Method for Simultaneously and Selectively Measuring AsIII, SbIII, and SeIV in Freshwater and Soils.

Anthropogenic and climatic perturbations redistribute arsenic (As), antimony (Sb), and selenium (Se) within the environment. The speciation characteristics of these elements determine their behavior and biogeochemical cycling, but these redox-sensitive species are challenging to capture, with few methods able to harmonize measurements across the whole plant-soil-ecosystem continuum. In this study, we developed a novel diffusive gradient in thin films (DGT) method based on aminopropyl and mercaptopropyl bi-functionalized mesoporous silica spheres (AMBS) to achieve in-situ, simultaneous, and selective quantification of AsIII, SbIII, and SeIV, three typical/toxic but difficult to measure inorganic species. When used for environmental monitoring within a river catchment, AMBS-DGT exhibited stable/accurate predictions of these species despite varying water chemistries (ionic strength 0.01-200 mmol L-1 NO3-, pH 5-9 for AsIII and SbIII, and pH 5-7.5 for SeIV). Furthermore, river deployments also showed that time-averaged species concentrations by AMBS-DGT were reproducible compared with high-frequency sampling and measurement by high performance liquid chromatography coupled with inductively coupled plasma mass spectroscopy. When AMBS-DGT was used for sub-mm scale chemical imaging of soil solute fluxes, the method resolved concomitant redox-constrained spatial patterns of AsIII, SbIII, and SeIV associated with root O2 penetration within anaerobic soil. Improved capabilities for measurement of compartment interfaces and microniche features are critical alongside the measurement of larger-scale hydrological processes that dictate the fine-scale effects, with the AMBS-DGT achieving this for AsIII, SbIII, and SeIV.

Validation and Assessment of Diffusive Gradients in Thin-Films (DGT) Technique for Measuring Nutrients in Taihu Lake Water with Algae Bloom.

DGT (diffusive gradients in thin films) technique has been developed for measuring nitrogen in freshwaters and applied to assess the bioavailability of phosphorus in soils/sediments. These two elements are the main nutrients causing algae bloom, but DGT has never been used in the field water conditions with algae bloom. In our study, a pair of DGT devices were used in comparison with grab sampling to characterize the performance of this technique to measure labile NO3-N, NH4-N, and PO4-P concentrations in algae-cultivated Taihu Lake water. The results showed that DGT measurement was highly affected by algae bloom and the environmental conditions using the current assemblies, especially for NH4-N measurement. For in situ measurement of nutrients in the real environment, an improvement to the DGT technique is required. The comprehensive assessment of the level of eutrophication needs to consider a variety of environmental factors rather than just the concentration of nutrients.

Molecular View on the Dissociation Pathways and Transactivation Regulation Mechanism of Nonsteroidal GR Ligands.

As a major drug target for anti-inflammatory therapy, the glucocorticoid receptor (GR) regulates a wide range of physiological processes through transactivation (TA) or transrepression. GR TA is involved in many adverse effects of GR-targeting drugs, and therefore, the discovery of novel GR ligands with lower TA activity and longer residence time is quite urgent. Undoubtedly, understanding the ligand dissociation mechanisms and the structural basis of the TA regulation is crucial for the development of novel GR-targeting drugs. Here, we used random accelerated molecular dynamics (RAMD) and funnel metadynamics (FM) simulations to explore the dissociation mechanisms of 5 classic glucocorticoids and 6 nonsteroidal GR ligands. Multiple ligand dissociation pathways were discovered. The classic glucocorticoids exhibit a strong preference for Path I, and most nonsteroidal ligands tend to dissociate along mixed pathways. We also find that the distinct unbinding preferences for AZD2906 and AZD9567, two representative nonsteroidal ligands with similar scaffolds but different TA activities, are primarily determined by their different polar interactions with the surrounding residues. Notably, the binding of AZD9567 poses a substantial impact on the conformation of the GR homodimer interface, which provides a valuable clue to understand the mechanisms of the TA-related side effects induced by the adjustments of the homodimerization process. These findings are critical for the structure-based rational design of novel GR ligands with more potent anti-inflammatory potency and reduced side effects.

Discovery of Novel GR Ligands toward Druggable GR Antagonist Conformations Identified by MD Simulations and Markov State Model Analysis.

Binding of different ligands to glucocorticoid receptor (GR) may induce different conformational changes and even trigger completely opposite biological functions. To understand the allosteric communication within the GR ligand binding domain, the folding pathway of helix 12 (H12) induced by the binding of the agonist dexamethasone (DEX), antagonist RU486, and modulator AZD9567 are explored by molecular dynamics simulations and Markov state model analysis. The ligands can regulate the volume of the activation function-2 through the residues Phe737 and Gln738. Without ligand or with agonist binding, H12 swings from inward to outward to visit different folding positions. However, the binding of RU486 or AZD9567 perturbs the structural state, and the passive antagonist state appears more stable. Structure-based virtual screening and in vitro bioassays are used to discover novel GR ligands that bias the conformation equilibria toward the passive antagonist state. HP-19 exhibits the best anti-inflammatory activity (IC50 = 0.041 ± 0.011 µm) in nuclear factor-kappa B signaling pathway, which is comparable to that of DEX. HP-19 also does not induce adverse effect-related transactivation functions of GR. The novel ligands discovered here may serve as promising starting points for the development of GR modulators.

Regulation Mechanism for the Binding between the SARS-CoV-2 Spike Protein and Host Angiotensin-Converting Enzyme II.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is mainly mediated through the interaction between the spike protein (S-pro) of the virus and the host angiotensin-converting enzyme II (ACE2). The attachment of heparan sulfate (HS) to S-pro is necessary for its binding to ACE2. In this study, the binding process of the receptor-binding domain (RBD) of S-pro to ACE2 was explored by enhanced sampling simulations. The free-energy landscape was characterized to elucidate the binding mechanism of S-pro to ACE2 with and without HS fragment DP4. We found that the stability of the T470-F490 loop and the hydrophobic interactions contributed from F486/Y489 in the T470-F490 loop of S-pro are quite crucial for the binding, which is enhanced by the presence of DP4. Our study provides valuable insights for rational drug design to prevent the invasion of SARS-CoV-2.

Pterostilbene inhibits hepatocellular carcinoma proliferation and HBV replication by targeting ribonucleotide reductase M2 protein.

Hepatocellular carcinoma (HCC), one of the most deadly diseases all around the world. HBV infection is a causative factor of HCC and closely associated with HCC development. Ribonucleotide reductase (RR) is a key enzyme for cellular DNA synthesis and RR small subunit M2 (RRM2) is highly upregulated in HCC with poor survival rates. We have previously shown that HBV can activate the expression of RRM2 and the activity of RR enzyme for the viral DNA replication in host liver cells. Thus, RRM2 may be an important therapeutic target for HCC and HBV-related HCC. Pterostilbene, a natural plant component, potently inhibited in vitro RR enzyme activity with the IC50 of about 0.62 µM through interacting with RRM2 protein, which was much higher than current RRM2 inhibitory drugs. Pterostilbine inhibited cell proliferation with an MTT IC50 of about 20-40 µM in various HCC cell lines, causing DNA synthesis inhibition, cell cycle arrest at S phase, and accordingly apoptosis. On the other hand, the compound significantly inhibited HBV DNA replication in HBV genome integrated and newly transfected HCC cells, and the EC50 for inhibiting HBV replication was significantly lower than the IC50 for inhibiting HCC proliferation. Notably, pterostilbene possessed a similar inhibitory activity in sorafenib and lamivudine resistant HCC cells. Moreover, the inhibitory effects of pterostilbine against HCC proliferation and HBV replication were significantly reversed by addition of dNTP precursors, suggesting that RR was the intracellular target of the compound. Finally, pterostilbine effectively inhibited HCC xenograft growth with a relatively low toxicity in nude mouse experiments. This study demonstrates that pterostilbene is a novel potent RR inhibitor by targeting RRM2. It can simultaneously inhibit HCC proliferation and HBV replication with a potential new use for treatment of HCC and HBV-related HCC.

The MCR-3 inside linker appears as a facilitator of colistin resistance.

An evolving family of mobile colistin resistance (MCR) enzymes is threatening public health. However, the molecular mechanism by which the MCR enzyme as a rare member of lipid A-phosphoethanolamine (PEA) transferases gains the ability to confer phenotypic colistin resistance remains enigmatic. Here, we report an unusual example that genetic duplication and amplification produce a functional variant (Ah762) of MCR-3 in certain Aeromonas species. The lipid A-binding cavity of Ah762 is functionally defined. Intriguingly, we locate a hinge linker of Ah762 (termed Linker 59) that determines the MCR. Genetic and biochemical characterization reveals that Linker 59 behaves as a facilitator to render inactive MCR variants to regain the ability of colistin resistance. Along with molecular dynamics (MD) simulation, isothermal titration calorimetry (ITC) suggests that this facilitator guarantees the formation of substrate phosphatidylethanolamine (PE)-accessible pocket within MCR-3-like enzymes. Therefore, our finding defines an MCR-3 inside facilitator for colistin resistance.

Assessing the potential availability of selenium in the soil-plant system with manure application using diffusive gradients in thin-films technique (DGT) and DOM-Se fractions extracted by selective extractions.

Knowledge of the Se fractionation and the role of dissolved organic matter (DOM) in soil is the key to understanding Se mobility and its bioavailability in the soil-plant system. In this study, single extractions using phosphate-buffer (PBS), sequential extraction procedures (SEP), and diffusive gradients in thin-films (DGT) were used to measure Se bioavailability in soil supplemented with selenite and organic amendment (cow and chicken manures). Selenium fraction was isolated into DOM-Se fractions, such as hydrophilic acid-bound Se (HY-Se), fulvic acid-bound Se (FA-Se), humic acid-bound Se (HA-Se), and hydrophobic organic neutral-bound Se (HON-Se), by a rapid batch technique using XAD-8 resin (AMBERLITE XAD™, USA). Simultaneous application of either cow or chicken manure with selenite could result in the decrease of Se availability in the soil. Isolating Se available fraction into DOM-Se fractions showed that low-molecular-weight DOM-Se as an available fraction and even HY-Se as a less available fraction (OM-Se) were likely the major sources for Brassica juncea (L.) Czern. et Coss uptake in soil. Moreover, knowledge of the DOM-Se composition, especially the low-molecular-weight DOM-Se fractions, is important for assessing the bioavailability of Se in soil, the results of which are more accurate than the chemical extraction method. The high value of Pearson correlation coefficients between CDGT-Se and Se concentrations in shoots, tubers and roots of Brassica juncea (L.) Czern. et Coss in cow and chicken manures treatment were 0.95 and 0.99, 0.96 and 0,96, and 0.89 and 0.97 (p < 0,05), respectively, indicating that DGT-Se can reflect the Se uptake ability by plants and can be used to predict the bioavailability of Se when manure and selenite are simultaneously applied.

The inflammatory cytokine IL-6 induces FRA1 deacetylation promoting colorectal cancer stem-like properties.

Colorectal cancer (CRC) has long been known for its tight association with chronic inflammation, thought to play a key role in tumor onset and malignant progression through the modulation of cancer stemness. However, the underlying molecular and cellular mechanisms are still largely elusive. Here we show that the IL-6/STAT3 inflammatory signaling axis induces the deacetylation of FRA1 at the Lys-116 residue located within its DNA-binding domain. The HDAC6 deacetylase underlies this key modification leading to the increase of FRA1 transcriptional activity, the subsequent transactivation of NANOG expression, and the acquisition of stem-like cellular features. As validated in a large (n = 123) CRC cohort, IL-6 secretion was invariably accompanied by increased FRA1 deacetylation at K116 and an overall increase in its protein levels, coincident with malignant progression and poor prognosis. Of note, combined treatment with the conventional cytotoxic drug 5-FU together with Tubastatin A, a HDAC6-specific inhibitor, resulted in a significant in vivo synergistic inhibitory effect on tumor growth through suppression of CRC stemness. Our results reveal a novel transcriptional and posttranslational regulatory cross-talk between inflammation and stemness signaling pathways that underlie self-renewal and maintenance of CRC stem cells and promote their malignant behavior. Combinatorial treatment aimed at the core regulatory mechanisms downstream of IL-6 may offer a novel promising approach for CRC treatment.

Reliability of Docking-Based Virtual Screening for GPCR Ligands with Homology Modeled Structures: A Case Study of the Angiotensin II Type I Receptor.

The number of solved G-protein-coupled receptor (GPCR) crystal structures has expanded rapidly, but most GPCR structures remain unsolved. Therefore, computational techniques, such as homology modeling, have been widely used to produce the theoretical structures of various GPCRs for structure-based drug design (SBDD). Due to the low sequence similarity shared by the transmembrane domains of GPCRs, accurate prediction of GPCR structures by homology modeling is quite challenging. In this study, angiotensin II type I receptor (AT1R) was taken as a typical case to assess the reliability of class A GPCR homology models for SBDD. Four homology models of angiotensin II type I receptor (AT1R) at the inactive state were built based on the crystal structures of CXCR4 chemokine receptor, CCR5 chemokine receptor, and δ-opioid receptor, and refined through molecular dynamics (MD) simulations and induced-fit docking, to allow for backbone and side-chain flexibility. Then, the quality of the homology models was assessed relative to the crystal structures in terms of two criteria commonly used in SBDD: prediction accuracy of ligand-binding poses and screening power of docking-based virtual screening. It was found that the crystal structures outperformed the homology models prior to any refinement in both assessments. MD simulations could generally improve the docking results for both the crystal structures and homology models. Moreover, the optimized homology model refined by MD simulations and induced-fit docking even shows a similar performance of the docking assessment to the crystal structures. Our results indicate that it is possible to establish a reliable class A GPCR homology model for SBDD through the refinement by integrating multiple molecular modeling techniques.